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Abstract 

Hybrid imaging techniques utilize couplings of physical modalities - they are called 
hybrid, because, typically, the excitation and measurement quantities belong to dif- 
ferent modalities. Recently there has been an enormous research interest in this area 
because these methods promise very high resolution. In this paper we give a review on 
hybrid tomography methods for electrical conductivity imaging. The reviewed imag- 
ing methods utilize couplings between electric, magnetic and ultrasound modalities. 
By this it is possible to perform high-resolution electrical impedance imaging and to 
overcome the low-resolution problem of electric impedance tomography. 

1 Introduction 

The spatially varying electrical conductivity, denoted by a — cr(x) in the following, provides 
important functional information for diagnostic imaging - it is an appropriate parameter 
for distinguishing between malignant and healthy tissue [31,38,48,49,52,65,75,91,98,113]. 
Although almost never stated explicitly, but very relevant for this work, the conductivity 
is frequency dependent, that is, a = c(x, w). The dependence on w is typically neglected, 
since the imaging devices operate at a fixed frequency. However, one should be aware that 
different image devices record conductivities in different frequency regimes. In figures la 
and lb, respectively, the specific values for a in healthy and malignant tissues are plotted, 
and it can be observed that cancerous tissue reveals a higher conductivity, in general. 




(a) Low frequency contrast of <r [38,52,98] (b) Radio-frequency contrast of a [48] 



Figure 1: Contrast of conductivity in biological tissue 

The standard approach for determining the electrical conductivity is Electrical Impedance 
Tomography (EIT). This approach is based on the electrostatic equation: 



V • (aVu) = , 



(1.1) 
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which describe the relation between the conductivity a and the electrostatic potential u — 
u(x). For imaging, EIT is realized by injecting series of currents at a certain frequency 
uj on the surface dfl of the probe and measuring the according series of voltages there. 
In mathematical terms the EIT problem accounts for computing a from the Dirichlet-to- 
Neumann mapping A CT : f i — y g, which relates the components of all possible pairs (/ = u\ gn , 
9 = (rd n u\ dn ). 

Precursors of EIT were for geological prospecting [45,95]. Current applications of EIT 
are in geology [2,85,89], but the applications in non-destructive testing for industry are of 
growing importance [106]. There has been considerable effort in adapting EIT for medical 
applications [18,23]. 

In mathematical terms, the EIT problem has been formulated first in 1980 by Calderon 
[21]. Uniqueness and stability have been studied extensively (see e.g. [1, 20, 101]). 1 The 
recovery of the electrical conductivity from the Dirichlet-to-Neumann map is heavily ill- 
conditioned, the modulus of continuity is at most logarithmic [3,74]. This result is consistent 
with the abilities of classical EIT to differentiate between conductors and insulators in 
geological prospecting (see table 1). On the other hand, the material differences in human 
and biological tissues are smaller, and logarithmic stability in this case guarantees only very 
low resolution. 



tissue type/material conductivity a 

skin (wet) 3 • lO" 3 S/m 

blood T-lO^S/m 

fat 2 • 10~ 2 S/m 

liver 5-10~ 2 S/m 

copper 6-10 r S/m 

granite (dry) 10~ 8 S/m 

Table 1: Conductivities of different type of tissue or materials, respectively, at 1 kHz [31], [28, Tab. 
9. If], [26] 



Recently, coupled Physics imaging methods (this is another name for hybrid imaging, 
which is more self explaining), have been developed, which can serve as alternatives for 
estimating the electrical conductivity. Passive coupled Physics electrical impedance imaging 
utilizes that electrical currents in the interior of a probe manifest themselves in different 
modality. In contrast, active coupled Physics imaging is by local excitement of the tissue 
sample. Mathematically, in most cases, coupled Physics imaging decouples into two inverse 
problems, which, in total, are expected to better well conditioned than the EIT problem [4]. 
Coupled Physics imaging led to a series of new results in inverse problems theory (see [13,56]). 

Coupled Physics imaging should not be confused with multi-modal imaging (actually both 
are sometimes called hybrid imaging), which is based on co-registering or fusing images of 
different modalities. In contrast to multi-modal imaging, the techniques described here, all 
involve interaction of physical processes, and therefore they have to use the correct modeling 
of the physics and chemistry involved. 

In this paper we survey mathematical models for coupled Physics electrical impedance 
imaging. We provide the principal mathematical modeling tools in section 2, which are 
various electromagnetic equations in different frequency ranges, basic equations of acoustics, 
and basics of magnetic resonance imaging (MRI). Based on this, we derive mathematical 
models for multiphysics imaging models in section 4. Bibliographic references are provided, 
covering experimental and mathematical research. 



The first step to prove unique recovery of a from Ao- for smooth conductivities, was the result in [99]. 
Now, there exist results for less regular conductivities (see the references in [f , f4.2.f .7] and [fOf , 6.f]). 
For planar domains, one even has uniqueness for L°°-conductivities [If]. 
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symbol 


name 


unit 




P(x) 


electrical charge density 


As/m 3 




E(x, t) 


electrical field strength 


V/m 




Hfx t) 


magnetic field strength 


A/m 




D(x,t) 


electrical displacement field 


As/m 2 




B(x,i) 


magnetic flux density 


Vs/m 2 = 


: T 


J(x,i) 


electrical current density 


A/m 2 




it(x,cj) 


electrical conductivity 


A/Vm = 




e(x,w) 


electrical permittivity 


As/Vm= 


=F/m 


ju(x,w) 


magnetic permeability 


Vs/Am= 


=H/m 




Table 2: Electromagnetic quantities 





2 Background 

In the following we provide basic background information on mathematical equations in 
Electromagnetism, Acoustics, and Magnetic Resonance Imaging (MRI). This overview is 
required to survey coupled Physics impedance imaging techniques in section 4. 



2.1 Electromagnetic theory 

We use as a starting point Maxwell's equations, which couple the physical quantities sum- 
marized in table 2. 



Maxwell's equations 

are the basic formulas of electromagnetism, and read as follows (see for instance [47] ) : 

dB 

VxE = -— , V D = p, 

VxH = J + — , V-B = 0. 

at 

For later reference, we recall two equations relating the fields E, B with the charge density 
p and the current density J: 

• First, the fields E and B act on free charges and electrical currents by the Lorentz 
force. The force per unit volume is given by (see [47, (5.12)]): 

f(x,£) = pE + J x B. (2.2) 

• Secondly, the work of the electromagnetic fields per unit volume and unit time is the 
power density [47, (6.103)], which in formulas expresses as 

iJ(x, i) = J • E. (2.3) 



In electromagnetism, in many applications, it is observed that the vector fields in (2.1) 
have the form 

E(x,i) = (£ fc (x)cos(u;t + ^!)) fc=li2) 3 > 
B(x,i) = (^(x)cos(wi + ^)) fc=123 , 
D(x,t) = (Z? fe (x)cos(wi + ^)) fc=12 3 , 
H(x,i) = (B*(x)cos(ujt + p k H )) k=12 3 , 
J(x,t) = (J fc (x)cosM + ^)) fc=123 . 
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Often the complex phasors 



E(x) = 


(# fc (x) e iuMpk z 


^ fe=l,2,3 


B(x) = 




) fc=l,2,3 


D(x) = 


(pl{x) e iuv » 


) fc=l,2 : 3 


H(x) = 


(# fc (x) e iu ^ 


) fc=l,2,3 


J(x) = 


(j fc (x) 


fe=l,2,3 



are used to rewrite (2.1) into a system of equations for complex functions without time- 
differentiation: 

VxE = -iwB, V D = p, 

. . . (2.4) 

V x H = J + IojU , V • B = 0. 

The real physical fields, which are time-dependent, can be uniquely recovered from the 
complex but time-independent phasors via the formula E(x, t) — Re (E(x) • e iLJt ). We do 
not distinguish between E(x, t) and E(x) in the following. 

For a variety of materials, including biological tissue, the vector fields E, D, B and H 
and J are connected by simple constitutive equations. In these linear media, the following 
material relations obtain: 

J = ctE . D = eE , B = /iH , (2.5) 

through which the material parameters a = er(x, t), e = e(x, t) and fi = /i(x, i) are defined. 
These are in general matrix-valued functions, depending on space and time. 

For biological tissue, the parameters conductivity and permittivity parameters have been 
investigated thoroughly [31,91]. These parameters are considered mostly as scalar, although 
rarely also anisotropy is investigated. A salient feature of the electrical parameters in tissue 
is their dispersion: That is, a = cr(x, uj) and e = e(x, uj) are functions of the frequency w 
of the electrical field, respectively [31,91]. However, as already stated above, this feature is 
mostly neglected, since the experiments usually operate in a fixed frequency regime. 

In biological tissue the electrical conductivity a varies from 1(T 4 to 10 2 S/m. The per- 
mittivity e and the permeability fj, are always scalar multiples of the fundamental constants 
e « 8,9-KT 12 F/m, ^ w 1,3-KT 6 H/m, respectively (see [47, Tab. A.3]). Moreover, in 
biological tissue, the ratio e/eo is between 10 and 10 8 [31]. Permeability is usually regarded 
as a constant [107, p. 151]. In the rest of the article, we assume that y, = (1q. 

For the introduction of potential functions, we consider the following two equations of 
Maxwell's system (2.4): 

V-B = 0, VxE =-iu)B. (2.6) 

Note, that we identified E, B and E, B, respectively. The first equation already guarantees 
that there exist a vector potential A satisfying 

B = V x A . (2.7) 

Consequently, by using (2.7) in the second equation of (2.6), it follows there also exists a 
scalar potential u, such that 

E = -iojA - Vu. (2.8) 

Through (2.8), the scalar potential u is determined by the observable field E up to an additive 
constant. Similarly, the vector potential A is not uniquely determined by the observable 
field B in (2.7). In fact, one is free to choose V • A (see [10, Thm. in 1.16]), e.g. 



V-A = 0. 



(2.9) 
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By choosing a particular constant for u and requiring (2.9), the two equations in (2.6) are 
equivalent to (2.7) and (2.8). - The requirement (2.9) is called Coulomb's gauge [47, 6.3], 
and we will use it in the rest of the article. 



2.2 Approximations for lower frequencies 

In the following we restrict attention to linear media, such as biological tissue. In mathe- 
matical terms this means that (2.5) are satisfied. 

Quasistatic approximation: 

We consider the third equation in (2.4), which is called the Ampere-Maxwell law. For linear 
media, it reads 

V x [ — B ) = crE + iuieE. 
VMo ) 

If uj satisfies 

w < a/e, (2.10) 

one can neglect the term iu>e~E (because it is small compared with crE). We replace it with 
the classical version of Ampere's law, 

V x ( — B ) = J. (2.11) 
\Mo J 

Using (2.11) together with the other three equations in (2.4), we get the following quasistatic 
approximation to Maxwell's equations: 

VxE = -!w/i(|H, V-(eE) = p, 

/ 1 \ (2.12) 
V x — B =ctE. V-B = 0. 

The quantities J := ctE and B appearing in (2.12) are related by the Biot-Savart law of 
Magnetostatics [47, 5.3]: 

B(x) = £2. / J(y) x j^^d 3 y, x e R 3 . (2.13) 
4?r J M 3 |x - y| J 

This law determines the magnetic field which is generated by a current density J. It is 
derived from the third and forth equation in (2.12), i.e. Ampere's law and V-B = (see 
for example [68, 5.1]). - Similarly, one can determine the vector potential corresponding to 
B (with Coulomb gauge) as 

f J(y) , 3 , 



A (x) = ^ / T^-rdty. (2.14) 
V ; 4ir J R3 \x - y| ( ' 

Certain dispersion models for a and s, like [31], suggest that the quasistatic condition 
(2.10) holds for a wide range of tissue types even at frequencies about 1 — 2 MHz. Note, 
however, that the available data on a and s have been obtained using in vitro measurements. 

At present, condition (2.10) is often also assumed because it simplifies the subsequent 
inverse problem of determining a. Whether the displacement current should be taken into 
the model or not depends on the specific application and the available data on a and e. For 
frequencies about 1 kHz, simulation studies have been conducted using a head resp. spine 
model [27, 102], in order to decide whether weE can be neglected. To the knowledge of 
the authors, there are at present no studies for higher frequencies, though further empirical 
and numerical research is necessary to determine the validity of the quasistatic regime in 
conductivity imaging. 
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Electrostatic approximation: 

We consider the first equation in (2.4), also called Faraday's law. In linear media, it reads 

V x E = -iuB. (2.15) 

It can be shown by a dimensional analysis (see [23, A. 2]) that if 

Vw7^CT<i?, (2.16) 

where R denotes a reference length of x, then Faraday's law can be replaced by the electro- 
static equation VxE = 0. From the latter equation, in turn, it follows that the electrical 
field can be represented as a gradient field, i.e., 

E=-Vu. (2.17) 

The electrostatic representation (2.17) is inserted into the complex time harmonic Ampere- 
Maxwell equation V x H = cjE + iweE from (2.4). Taking the divergence, one gets 

V-(kVu) = 0. (2.18) 

Here, the quantity 

«;(x, a;) = <r(x, to) + icj e(x, oS) (2-19) 
is called complex conductivity and is measured in S/m as well. 

Often, though, one combines the quasistatic with the electrostatic approximation: If we 
apply the divergence to J in the quasistatic version of Ampere's law (2.11), this yields the 
real- valued elliptical equation in (1.1), 

V • J ( = 5) -V • <tE (2 - ?) -V • (o-Vm) = 0. (2.20) 

We use later that in the quasistatic case, the magnetic field generated by the current J = crE 
is given by the law of Biot-Savart (2.13). Similarly, the vector potential generated by J can 
be calculated by (2.14). 

Experiments suggest that for u> larger than 1 MHz, the electrostatic simplifications are 
not valid anymore and it is reasonable to use a more accurate model like the quasistatic 
approximation or the whole set of Maxwell's equations [96] . 



2.3 A reduced eddy current model 

In this section we derive a simplified mathematical model of the full Maxwell equations (2.4) 
resp. the quasistatic Maxwell equations (2.12). Such a model was developed in [34] to de- 
scribe eddy currents in biological tissues (see also [29]), with the prime goal for improvement 
of induced current electrical impedance tomography (ICEIT) . 

Here, we assume knowledge of a current distribution with compact support in a set 
F, where F is a bounded domain outside the conducting material in the domain Q. For 
instance in ICEIT, is the excitation current, and F is the domain of the excitation coil. 

We start from the time-harmonic Maxwell equations (2.4) in the source-free case p = 0. 
These are simplified below. With an imposed current density Ji, the equations read as 
follows 

VxE = -i^ H, V-(eE) = 0, 

/ l \ (2.21) 
V X I — B ] = kE + V-B = 0, 

VMo / 

where we use the complex conductivity k as in (2.19). 
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We separate the electrical field E into primary and secondary components, and describe 
those via their potentials (cf. (2.8)): 

E n = — iujAr, 

E = E D + E s with P P (2.22) 

p E s = -iujA s -Vu. V ' 

Here, A p is the vector potential generated by the imposed current (see (2.14)). A p is 
assumed not to depend on the tissue conductivity. The secondary field A s is the vector 
potential generated by crE. 

For the representation of E p we do not need a scalar potential, whereas for E s we do. 
Indeed, assume that E p = — iujA p — Vu p . Because the permittivity e is constant in vacuum, 
we have V • (eE p ) = eV • E p = (second equation in (2.21)). We use Coulomb's gauge (2.9) 
for A p , yielding 

= V • E p = V • (-i(jA p - Vu p ) = -iu)V ■ A p - V • Vu p ( = 9) -Au p . (2.23) 

In the source-free case p — 0, the only solutions making physical sense are constant, and 
without loss of generality, we can take u p = [47, 5.4,6.3]. In particular this means that 
Ep is completely determined by the vector potential A p . 

Now for biological tissue, it was observed in [34] that kE w — K(iujA p + Vu). Therefore 
we can work with 

E"™ := -iujA p - Vu (2.24) 

as an approximation for the electrical field. Using the third equation in (2.21), we have 
V • (V x B) = V ■ J = 0. Consequently, the following condition holds for the scalar potential 
u: 

V • (kVu) = -iujA p ■ Vk (2.25) 

At the boundary dfl of the conducting region we have ftE™™ • n = 0, as there is no current 
flux through the boundary. This yields the condition 

d n u = —iujA p - (2.26) 

at the boundary. Note that in the later applications, the equation (2.3) will only be con- 
sidered on a bounded domain Q. To ensure uniqueness for the problem (2.3), (2.26), one 
usually requires a scaling condition for u, e.g. J gn u dS — 0. 

With Ap (calculated from the current Jj) and it, we find the electrical field ~E new and 
a magnetic field B™ et0 (from (2.14)), which approximately satisfy (2.21). This has been 
checked in numerical simulations [34, A.C], [29] (see also [25]). 

For practical simulation and reconstructions, one often restricts oneself to the quasistatic 
version of the eddy current model. That is, one neglects the imaginary part in k and can 
therefore work with the real electrical conductivity a instead. Then one replaces with the 
model: 

V • (aVu) = -iujA p ■ Ver (2.27) 

As noted in the discussion of the quasistatic condition (2.10), which model is to be used de- 
pends on the specific application and the available empirical data on the material parameters 
a and e. 



2.4 Acoustics 

Equations for the acoustic pressure are obtained by linearization from the fundamental 
equations of fluid dynamics, which relate the pressure p — p(x, t), the density p — p(x, t), 
and the velocity v = v(x, t) of the fluid, respectively. 

Because it is a linearized model the variations of these parameters relative to a ground 
state (po( x ), Po( x ), 0) [64, §64] have to be small to be consistent with reality. 
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In particular, the standard wave equation 



□ p:= \dttP- Ap = 0, (2.28) 



c 2 



is derived from the linearized conservation principles for impulses, mass, and a relation 
between pressure and density, which read as follows: 

p a d t v + Vp = 0, (2.29) 

d t p + p V-v = 0, (2.30) 

p = c 2 p. (2.31) 



Here, c = c(x) = y / d p p(po(xj) denotes the speed of sound in fluid dynamics (see [64, § 64]). 
Applying the divergence to (2.29) and calculating the time derivative of (2.30), we eliminate 
v to get Ap — duP- Substituting the time derivative of (2.31), we obtain the wave equation 
(2.28). 

In the following, we assume that the fluid is perturbed by electromagnetic fields, causing 
the Lorentz force (2.2) and Joule heating (2.3). We incorporate these effects into the equa- 
tions for the ultrasound pressure p. We first observe that the Lorentz force density f enters 
as a source term into the balance equation for the impulse (2.29), and we get 

Po d t v + Vp - f = 0. (2.32) 

Secondly, the linearized expansion equation, describes the relation between the energy ab- 
sorption, described by the power density function H, and the pressure 

d t p = YH + c 2 d t p . (2.33) 

Here Y is the dimensionless Griineisen parameter. Very frequently it is assumed constant, but 
in fact is material dependent, and thus Y = Y(x). (2.33) can be derived from thermodynamic 
relations and the principle of energy conservation [40, 0.1]. We emphasize that (2.33) is a 
generalization of (2.31). 

The wave equation with outer force and energy absorption is now derived from (2.32), 
(2.33) and the mass conservation principle (2.30). We arrive at the following inhomogeneous 
wave equation 

□ p = -V-f(x,t) + ^a t if(x,f) (2.34) 

We will use this equation with f(x, t) — S(t)(f>(x) and H(x,t) = S(t)ip(x) in sections 4.2.1 
and 4.2.3. 

Remark 1. By DuhameVs principle (see [30, p. 81]), the solution of 

n P = 8(t)(t>(x)+d t 5(t)il>(x), (2.35) 

where S is the S- distribution, is the sum of the solutions of the inhomogeneous wave equations 

□ Pi(x,t) = □P 2 (x,t) = 

Pi(x,0) - P>(x,0) = V(x) (2.36) 

atPi(x,0) = 0(x), 5 t P 2 (x,0) = 0, 



respectively. 
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3 MRI imaging 

General comments and modeling of MRI 

Magnetic Resonance Imaging (MRI) is frequently used as the basis for coupled Physics 
imaging. 

Since the modeling of MRI is not so common in Inverse Problems, it is reviewed here: 
MRI visualizes the nuclear magnetization M = M(x, t) resulting from selectively induced 
magnetic fields B = B(x, t). These two quantities are related by the simplified Block equa- 
tion 2 

d t M. — 7M x B. (3.1) 

Here 7 denotes the gyromagnetic ratio specific to a proton. For instance, for hydrogen one 
has 7 = 4.6 • 10 7 T"^- 1 . (T ... Tesla) 

For a stationary magnetic field in z-axis direction, B(x, t) = B(x) = (0,0, B(xj), the 
solution of the simplified Bloch equation with initial data M(x, 0) = Mo(x) is given by [37, 
2.3.2] 

/ M§ cos(-7-B(x)i) - Ml sm(-jB(x)t)\ 
M(x, t) = [ M§ sin(-7.B(x)t) + M% cos(- 7j B(x)t) . (3.2) 
V MS J 

This identity shows the rotation of the magnetization M around the z-axis when applying 
the stationary magnetic field (0, 0, B(x)). The quantity 

w(x) = -7-B(x) (3.3) 

is called the Larmor frequency, and is an essential ingredient for MRI. 

We explicitly introduce two particular magnetic fields and their Larmor frequency: 

• For the static field 

B := (0,0, B ), (3.4) 
the Larmor frequency is denoted by 

loq := -7-80- (3.5) 
For the case of hydrogen, and with |£>o| = 1.5 T, one has |wo| = 63.9 MHz. 

• By a gradient field corresponding to the vector G we mean the field 

B G (x) := (0,0, G-x). (3.6) 
If B G is applied along with Bo, the Larmor frequency is w(x) = ojo — 7G • x. 

Let now B(x, t) — B(x) be a static magnetic field and M(x, t) the resulting magneti- 
zation due to (3.1). The MRI signal is collected from voltage measurements of a detector 
coil. 

It is known that the voltage V(t) received in a detector coil satisfies 



V(t)oc J M(x,<)-B roc (x)d 3 x, 



where B rcc is the Biot-Savart-field (2.13) corresponding to a unit current in the detection 
coil [37, 7.2]. 



2 We restrict our discussion to the influence of B on the precessional motion of M. More adequately, the 
motion of M depends on the material-specific relaxation paramters Ti, T2 [37, 4.4], For the discussion 
of hybrid techniques, it suffices to study the simplified Bloch equation, and subsequently, the simplified 
MRI signal in (3.7). 
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Bo 



B g 2 
B G - 



Figure 3: Time shape of a pulse sequence with three different gradient fields B Gx , B g h, B g *. The signal 
is measured during application of B Gx . 

Let us assume that B rcc = (r cos(#b), r sin(0s), 0). Then using (3.2) and trigonometric 
identities, one finds that the voltage is 



V(t) =K J M sin(u;(x)f + 6 B - Md 3 x. 



Here, Mrj(x) and (f>o{x) are the polar coordinates of the transverse magnetization in the 
complex representation 

Af e**° = M$ + iMg. 

Here, the terms on the right hand side are the x- and y-components of the magnetization 
vector M = (Mg , M 9 , M Z ) . 

In a signal analyzing step, one multiplies this voltage with the reference fields sin(wo t) 
and cos(wo t) respectively. After filtering, ond finds the real and complex part of the following 
quantity (up to a fixed complex proportionality factor K): 

S(t) = I ^(xJe^We^W'At (3.7) 



S(t) will be referred to as the MRI signal [37, 23.7]. 

For MRI imaging, we need only the expression for the signal for the case of a a static 
field B(x) as in (3.7), because the signal is only measured while a static field is applied (see 
figure 3). 

The signal S(t) depends both on w(x) as well as on M. In actual experiments in MRI, 
both quantities arc varied to collect data. 

Experiment: 

In magnetic resonance imaging experiments one applies different fields B in order to vary 
w(x) and M, such that sufficient data of S(t) can be collected to image the interior of the 
specimen. In a general pulse sequence of MRI is a superposition of three magnetic fields, 
Bo, Bi, B G , 

B(x, t) = B (x) + Bi(x, t) + B G (x, t) . (3.8) 

B and B G are as in (3.5) and (3.6), respectively. Moreover, Bi = (Bi(t) cos(cJii),.Bi(i) sin(wii 
is a radio- frequency field. The time-shapes of the applied fields are depicted in figure 3. 

The field Bi is responsible for creating a nontrivial Mq as initial data: At every point 
where the resonance condition 

uji = w(x) (3.9) 
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holds, the magnetization vector is flipped away from the z-axis about a certain angle, and 
inconsequence, a non-zero transverse magnetization Moe 1 ^ results. 

Gradient fields with different vectors G are induced at different time (see figure 3). 
Applying B G ~ with G z = (0, 0, G z ) during the radio-frequency field, the resonance condition 
(3.9) is restricted to a small slice around the plane z = (and only there Mq ^ 0). Every 
signal (3.7) recorded subsequently can be approximated by a 2D integral (cf. [37, (10.34)]). 

One applies B Gy with G y = (0, G y ,0) between the radio-frequency pulse and the mea- 
surement for a certain time t' . At this time, the Larmor frequency changes to w(x) = 

+ G y y. Afterwards, the phase of the transverse magnetization has changed from cf>o to 
0o + t'G y y. 

Now after making these preparations, measurements are taken and the signal (3.7) is 
recorded. During the measurements, though, a gradient field B Gz with G x = (G x ,0,0) is 
applied, which changes the Larmor frequency in (3.7) to w(x) = ujq + G x x. 

The upshot of applying these three types of gradient fields is that one can modify 4>o and 
w(x) in (3.7), such that one has access to the data 



Introducing the 2D vector k = (tG x ,t'G y ), (3.10) is just the 2D-Fourier transform with 
respect to k of the transverse magnetization Mo(x, y, 0)e^° , sampled at the frequency vector 
k. By repeating the pulse sequence and varying (G x ,G y ), and thus k, one collects data in 
the frequency regime from the Fourier transform of Mq(x, y, 0)e 4 ^°, 



The Fourier transform allows to recover Mo(x, y, 0)e^° . 

Remark 2. The modeling assumes that the magnetic fields in (3.8) have an ideal shape, 
and are not affected by the body tissue, in particular. 

In fact, the radio- frequency field Hi is affected by electrical properties of the tissue. This 
effect can be utilized to reconstruct the conductivity cr(x, o>o) using the Larmor frequency 
ujq [50, 51, 111]. - This is one of several examples of imaging the electrical conductivity, or 
the current density at the Larmor frequency, using the measurement setup of MRI (for other 
approaches see [83, 103]). 

4 Coupled Physics models for conductivity imaging 

In this section we review coupled Physics techniques for imaging the scalar conductivity a 
from (2.5). 

4.1 Magnetic coupling 

Here, we discuss two sorts of methods for conductivity imaging based on magnetic cou- 
pling. One set of methods uses coupling of magnetic and electrical phenomena described 
by Maxwell's equations. The other set of methods uses MRI and exploits the influence of 
electrical currents on the MRI signal. Both kinds of methods have either contactless exci- 
tation or contactless measurement (or both). The latter is a practical advantage in clinical 
experiments. 

4.1.1 Magnetic resonance-EIT (MREIT) / Current density impedance imaging 



In MREIT and CDII, respectively, one perturbs the MRI signals by injecting an electrical 




(3.10) 




(3.11) 



(CDII): 



current. 
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Experiment: 

The measurement setup consists of an MRI machine and, additionally, electrodes on the 
surface of the specimen (within the MR machine). The MRI pulse sequences are performed, 
and an electrical current is injected for a short period, before the induction measurement. 
The injection perturbs the Larmor frequency and leads to a change in the nuclear magneti- 
zation, which alters the MRI signal. In CDII, in addition, the experiment is repeated with 
the probe rotated relative to the MRI machine. 



Mathematical modeling: 

The mathematical modeling of these two imaging techniques is very similar to the modeling 
of MRI in section 3. The difference is that the injected electrical currents manifest themselves 
in an additional magnetic held, namely, the Biot-Savart held Bj = (Bj , Bj, Bj) from (2.13). 
This field has to be taken into account in (3.8), in the mathematical model of MRI. Thus 
the total magnetic field in MREIT/CDII is given by 



B(x,f) = B + Bi(x,i) + B°(x,t) + Bj(x,t) . 



The current J is injected before the induction measurement. As J is a direct current or has 
very low frequency, the Biot-Savart field Bj acts like an additional gradient field B * . As 
discussed in section 3, this results in a change of the Larmor frequency Aw = — 7|Bj|, at the 
time Tj when the current J is applied. The transverse magnetization therefore accumulates 
an additional phase difference A(f> = 7|Bj|Xj, which in this situation is approximately 
jB^Tj (see [37, 27.3.6]). 

Compared with the signals (3.11) recorded in MRI, the signals in MREIT/CDII reveal a 
phase shift. Because in the experiments two currents of opposite sign are applied [93,94,108], 
one collects the following data: 

S±(k x ,k y )= f M o (x,y,0)e^° e ±*r*5( a ,i/,o)7b e * k ' x d(x, y), (4.1) 

JR 2 

Here, ± indicates two currents of opposite signs. Using the inverse Fourier transforms of 
(4.1), and denoting 

m± = M(x, y, zq) e ±l ~< B ^ x > y ^ Tj , (4.2) 



one finds that: 



^(x ;y ,,o)=^i m (iog(^±^4n 

2 7 Tj V \m-{x,y,z )J 



which is the component of Bj in direction of the main magnetic field Bo- 

Other components of Bj can be obtained by rotating the sample relative to the MRI 
machine. Then one determines the current density by Ampere's law (2.11) applied to B = 
Bj, that is, J = x Bj. 



Imaging: 

B J and J depend non-linearly on the conductivity a. Depending on the acquired MRI data 
it is common to differ between MREIT, that uses values of B z , and CDII, where knowledge 
of the current density J is assumed. CDII, in order to obtain the three dimensional vector 
B J , requires rotating the probe relative to the MRI-machine three times. 

Problem 1 (CDII). For a given domain O, representing the domain of the probe to be 
imaged, let boundary data f, g be given. Moreover assume that on the whole of ft, — V x 
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Bj = J 



trVu, or |J|, respectively, are given, where u and a are related by 

V • (o-Vu) = in n C K 3 
u = f on <9f2 
a d a u = g on d£l. 



(4.4) 



The imaging problem consists in determining a. 

Problem 2 (MREIT). Let boundary data f , g be given. Let a and u be related by (4-4)- 
The imaging problem consists in determining a from B j , which is related to a and u via the 
Biot-Savart relation (2.13), 



MREIT/CDII have been suggested in the late 1980ies [92]. Since then, efficient algo- 
rithms have been developed for solving the two problems. There also exist a variety of 
experimental biomedical studies [94,108]. In practice, the experiments are repeated at least 
two times with different boundary data, which provides more independent information on 
a. 

Mathematics of MREIT and CDII: 

For CDII, many uniqueness and stability results have been obtained: 

• Let |Ji|, | Ja ] and corresponding Neumann data a d n u = gi, i = 1,2 be given, 
where the currents satisfy Ji Jjf J2 everywhere in fi. Then, a can be recovered 
uniquely [53,54,63]. 

• For two given vector-valued current densities, there exists an analytical recon- 
struction formula, which determines a up to a multiplicative constant [41,67]. 

• With one current density vector J and Dirichlet data u\gn given, unique recovery 
of a for the two-dimensional version of problem 1 is guaranteed [62] . 

• Recently, it was observed that the absolute value of one current, |J|, and Dirichlet 
data [79,80] are sufficient for recovering a. These uniqueness results for CDII have 
been proven by interpreting level sets of u as minimal surfaces or by exploiting an 
energy-minizing property of u [78-81]. They require u to be free of singularities. 
Stability results of the |J| based reconstruction have been given by [78,82]. 

For MREIT: 

• Uniqueness and stability results for MREIT are given in special cases by [71,72]. 

• Various numerical algorithms have been developed, which reconstruct a using 
Bj. and Neumann data <?j, i > 2. For an overview, see [94, Sec. 3]. 

4.1.2 Magnetic detection impedance tomography (MDEIT): 

MDEIT uses magnetic boundary detection to recover currents and conductivity inside the 
tissue. 

Experiment: 

The equipment are electrodes to inject currents and several coils to record the magnetic 
field. A current is injected into the probe and the strength of the magnetic field is measured 
by detector coils at different positions of the boundary. This technique is discussed in [46,55] 
and for practical applications see [73]. 
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Mathematical modeling: 

The mathematical model is the electrostatic model, V • (erVu) = from (2.20), which relates 
the electrical potential u and the conductivity a. Neumann boundary conditions are given by 
the injected currents. The Biot-Savart law (2.13) provides the relation between J = —a'Vu 
and the magnetic field B measured by the detector coils. 

Imaging: 

The imaging problem can be decomposed into two subproblems. In the first part, one 
reconstructs the current density J inside $1 from measurements of B at the boundary dQ. 
In the second part, one faces exactly the same problem as in CDII, namely to reconstruct a 
from interior data of J (problem 2). 

The inverse problem in the first step is a linear one: to recover J given the Biot-Savart 
law. We can restrict ourselves to current densities satisfying V ■ J = 0, see (2.20). 

Problem 3 (Inverse problem MDEIT, first step). Determine J|n from B|an using the Biot- 
Savart relation 

B(x) = ^ / J(y) x ^^U 3 y, x G R 3 , (4.5) 

taking into account that V • J = 0. 

The second problem is analogous to CDII (problem 2): 

Problem 4 (Inverse problem MDEIT, second step). Let be given data f and g on dfl. 
Determine a from J = — ctVu, which are related by (4-4) ■ 

Mathematical analysis: 

The mathematical investigations relating to the linear problem 3 concern uniqueness is- 
sues and development of numerical algorithms. The kernel of the Biot-Savart operator, has 
been investigated in [42,55]. Uniqueness in the reconstruction problem has been shown for 
particular cases, such as directed currents in layered media [42,55,87]. Different regular- 
ization methods have been applied to solve problem 3, and tested using real and simulated 
data [43,46]. 

4.1.3 EIT with induced currents (ICEIT) 

ICEIT excites currents by magnetic induction and generates a voltage, which is measurable 
at the boundary of the probe. 

Experiment: 

For the experiments, one uses coils and electrodes. A time harmonic current is produced in 
the excitation coil and eddy currents are induced in the tissue. The voltages are detected 
at the boundary of the probe [34,112]. 

Modeling: 

We recall the setting of the reduced eddy current model in section 2.3, referring to figure 
2. The excitation coil lies in a domain F where the excitation current J,; is supported. The 
conducting tissue is in the domain f2, disjoint from F. 

In (2.24), we have introduced a reduced version of the electrical field as 

E = — iuA p — Vtt, 
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where A p is the primary vector potential generated by the excitation current (see (2.14)). 
The scalar potential is governed by (2.27) with boundary condition E • n = on dil. To 
ensure a unique solution, we require that J gn udS = 0. This gives the description of the 
eddy currents inside the tissue. 

The measurement data consist of (finitely many) line integrals E • dl. Because the 
determinatino of A p via (2.14) does not depend on inhomogeneities of a, one obtains line 
integrals of Vtt, which are potential differences. As in classical EIT, these correspond 
approximately to knowledge of Dirichlet data u\qq, if one chooses a voltage reference. For 
this, we assume J gQ u dS = 0. 

Imaging: 

The imaging problem in ICEIT is an inverse boundary value problem: 
Problem 5 (Inverse problem ICEIT). Let (<r,u) satisfy 

V • (<rVu) = — iwAp • Vcr in 

(4.6) 



d a u = — iwAp • n on dQ 



udS = 0. 

an 



Given are boundary data u\gn and the primary potential A p |q in Q, which is related to the 
excitation current Ji on M 2 by 

Mo f Ji(y) ,3, 



Determine a . 



A P (x) = ^ / P^TdV (4.7) 
4tt J R 2 |x - y| 



Mathematics and Numerics: 

ICEIT has been investigated numerically: Different regularization procedures have been 
applied, e.g., in [34,112]. 

4.1.4 Magnetic induction tomography (MIT): 

In MIT, eddy currents are induced, and the resulting magnetic fields are measured at the 
boundary. 

Experiment: 

The measurement equipment consists of coils placed at the boundary of the conducting 
tissue. Electrical currents in the coils induces eddy currents in the tissue. These currents 
generate a secondary magnetic field, which is measured by detection coils [36,97], outside of 
the probe. 

Modelling: 

The mathematical model is based on the reduced eddy current model (2.21) in section 2.3, 
as in the previous section. 

By (2.24), the electrical field is represented as E = — icjA p — Vit, where A p is the primary 
vector potential, generated by the excitation current Ji (see (2.14)). The scalar potential u 
is determined by (2.27). 

The signal recorded in a detection coil is the electrical voltage induced by E|n- Special 
types of detectors can be constructed, which are not influenced by the primary excitation 
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field E p = —icuAp, but only by the reaction field [36]. As in (2.22), the electrical reaction 
field is E s = — icjA s — Vw, where A s is the vector potential generated by the current crE 
(see also [29, (25)]). 

The measured voltage in a coil with surface C is then 

v = f E s • dl = / n • (V x E s ) = -iu / n • B s . 

Jac Jc Jc 

In the last equation, we have applied Faraday's law (2.15) to E = E s and B = B s . B s is 
the magnetic field generated by the current current 

<rE|n = — icrcoAp — a'Vu , 

see (2.13). 



Imaging: 

The reconstruction problem of MIT is as follows: 

Problem 6 (Inverse problem MIT). Let A p be the vector potential generated by the exci- 
tation current Ji, given by (4.7). Let u satisfy the Poisson equation (4.6) with parameter 
a. 

The given data are different voltage measurements 

V k = f n ■ B s 

Jc k 

acquired with detector coils on the surface Ck- Here, B s is function of a, u and A p as 

B s (x) = (<wA p (y) + V«(y)) x * ~ ^ 3 d 3 y, x e R 3 . (4.8) 

4vr J Q \x-yr 

The imaging problem consists in determining a . 



Eddy current imaging techniques play an important role for non-destructive evaluation in 
industry [100]. Actually, magnetic induction tomography has been developed for industrial 
non destructive evaluation and geological surveying [19,77,86]. For these application the 
reduced model of magnetic induction is not valid. However, also biomedical applications 
of MIT have been proposed later and the reduced model is now in general use. But MIT 
does not provide high resolution [36]. The primary importance of biomedical MIT is for 
contactless monitoring of physiological processes such as respiration [97]. However, it was 
observed recently that coupling MIT with ICEIT (sect. 4.1.3 and classical EIT leads to 
better resolution of the corresponding linearized inverse problem [35] . 

4.1.5 MREIT with magnetic induction (ICMREIT) / Current density impedance 
imaging with magnetic induction (CDII-MI): 

In ICMREIT and CDII-MI, one perturbs the MRI signals by an electrical current, which is 
produced by magnetic induction. 

Experiment: 

The measurement apparatus is an MRI machine with one or more coils for inducing eddy 
currents in the probe. MRI pulse sequences are performed. Additionally, eddy currents 
are produced before the induction measurement, and the perturbation of the MRI signal is 
recorded. 

The technique has been first proposed in [84]. 



HYBRID TOMOGRAPHY FOR CONDUCTIVITY IMAGING 



17 



Modeling: 

The modeling of these techniques is similar to MREIT/CDII discussed in section 4.1.1, 
except that we use the reduced eddy current model (2.27) from section 2.3 to model the 
electric field. 

By (2.24), the eddy current inside the tissue is J — —a(iujA p + 'Vu), with A p determined 
as vector potential generated by the impressed current density . This (low frequency) eddy 
current generates a magnetic field component Bj as in (4.8). 

As discussed in the section 4.1.1 on MREIT/CDII, the signal measured has information 
of one component of the magnetic field, as in (4.3); this information is used in ICMREIT. 
By rotating the probe, one can obtain the full vector field Bj, which is used in CDII-MI. 

Problem 7 (Inverse problem CDII-MI). Let A p be the vector potential generated by the 
excitation current Ji, given by (4.7). Let u satisfy the Poisson equation (4.6) with parameter 
a. 

Given are data J = — cr(iwA p + Vu). Determine a. 

Problem 8 (Inverse problem ICMREIT). Let A p be the vector potential generated by the 
excitation current Ji, given by (4.7). Let u satisfy the Poisson equation (4.6) with parameter 
a . 

Given are data Bj with Bj, as in (4.8). Determine a. 
Mathematics and Numerics: 

Algorithms for MREIT/CDII have been adapted for the numerical solution of these prob- 
lems. Thereby using measurements involving Bj [33,84], or measurements involving J [72]. 

4.2 Coupling involving acoustic modalities 

Up to now, we have discussed how the conductivity distribution behaves under influence of 
electric and magnetic fields. In the last ten years, there have been several proposals to utilize 
couplings between electrical and acoustical properties. Ultrasound propagation translates 
well from outer measurements to inner quantities of the probe. It has long been used well 
for other hybrid techniques such as photo-acoustic imaging. Therefore high resolution can 
be expected. We explain the proposals in the following, and draw especially on background 
in section 2.4. 

4.2.1 Impedance-acoustic tomography (IAT): 

Impedance-acoustic tomography exploits Joule heating due to electric currents to produce 
an ultrasound wave. 

Experiments: 

The measurement equipment are electrodes and ultrasound transducers. The electrodes 
are attached to the surface of the tissue, and an electric current is injected into the tissue. 
Joule heating results, and the tissue emits an ultrasound wave due to thermal expansion. 
Ultrasound transducers record this pressure wave at the surface of the probe or in a contact 
medium, e.g. water. The method has been proposed in [32]. 

Modeling: 

We use the quantitative description of the ultrasound wave incorporating electric effects 
from section 2.4. Since there is no additional outer force f , the relevant wave equation is 
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iflttpfot) - Ap(x,t) = Jdj-H", 

The model has to be supplied with the initial data p(x, t) = and ■§^p(x,t) = 0. The 
equation can be considered in 1" x (0, oo), if one notes that H has its support in tt. 

We assume that H(x.,t) = H(x)g(t), and that g(t) can be assumed as delta pulse. This 
is because the time of the electrical pulse multiplied with the sound speed in the tissue is 
negligible. We therefore obtain an inhomogeneous wave equation as outlined in remark 1 in 
section 2.4: 

-^d tt p-Ap=0, (x,f)e K n x(0,oo), 

p =A ct |e| 2 , (x,i)e R"x{0}, ( 49 ) 

Up 

d tP = 0, (x,i)e R"x{0}, 

where a vanishes outside the domain CI. 

The power density H can be related to the conductivity with Joule's law (2.3), that is, 

H = J-E = cr|E| 2 , (4.10) 

with unit W/m 3 , quantifying the energy converted per unit volume and unit time. The 
power density H is the intermediate link between the electrical excitation and the acoustical 
detection, and it appears in the wave equation as a source term. 



Imaging: 

Recovery of the initial pressure from boundary data is a long-established technique. Re- 
construction is possible with different measurement geometries, and it has been mainly 
developed for photo-acoustic tomography, a hybrid technique combining laser excitation 
and ultrasound detection [59,69]. 

If we assume that the quantities j3 and C p are constants, then this gives knowledge of 
the power density er|E| 2 . 

The second step is the reconstruction of conductivity from the interior dataer|E| 2 = J E 
throughout il. Obviously, we have to indicate how the electrical field entering this problem 
as initial source state is to be described. Here it makes a difference which kind of frequency 
is employed. 

In the electrostatic case of impedance-acoustic tomography [32] , we can use the approx- 
imation (1.1) known from classical impedance tomography, including Dirichlet data. Then 
we end up with the following inverse problem with interior data: 

Problem 9 (IAT, conductivity reconstruction). Let 

V • (aVu) = in O , 
u\dn = f ondQ, 
a d n u = g on dfl . 

Given are f, g, as well as interior data H = cr|Vu| 2 . Determine the conductivity a from 
these. 

Mathematical analysis: 

Unique recovery of a was shown in [22] for the two-dimensional version of problem 9, us- 
ing Hxx, H\2, H22 for Hij — ctVwi • Vit2- (The term H12 can be obtained using ct|Vui| 2 , 
er|Vui| 2 , <t|V(ui + «2)| 2 by polarization.) This result was generalized in [14] to three 
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dimensions, using multiple measurement data to guarantee unique global reconstruction. 
Reconstruction of a using one power density H = er|Vu| 2 , based on the nonlinear PDE 

V • (H .J?" 1 = was investigated in [121. Numerical reconstruction procedures have been 

| Vu| J 

applied to simulated data in [22,32,58]. The stability result from [60] applies to reconstruc- 
tions of a from two or more power densities. 

4.2.2 Quantitative thermo-acoustic tomography (qTAT): 

In thermo-acoustic tomography, one exploits Joule heating due to microwave excitation to 
produce an ultrasound wave. 

Experiment: 

The measurement equipment consists of a microwave generator and ultrasound transducers. 
Microwaves pulses are generated with a frequency in the GHz range. They propagate within 
the tissue sample and energy absorption causes Joule heating. The tissue reacts with a short 
expansion wave which is recorded outside the tissue, either on its surface or in a contact 
medium such as water. Measurement setups are described, e.g., in [69]. 

Modeling: 

The general principles of the technique are the same as for impedance- acoustic tomography, 
described in section 4.2.1. We use again the inhomogeneous wave equation (2.36) with 
f = 0. Assuming that the power density is of the form H(x, t) = H(x)g(t) with g(t) = S(t) 
a delta pulse. Therefore, we arrive at the same wave equation (4.9) with initial pressure 
F(x)=J(x).E(x)= ( x(x)|E(x)| 2 . 

A high-frequency description is necessary for modelling the oscillating field E evoked by 
microwaves. Here one starts from the time derivative of the Ampere-Maxwell law in (2.1), 

V x d t B = Li d t J + Mo<9 t fD - S, 

and couples it with Faraday's law 9(B = — V x E to get the radiation equation 

e^o^ttE + <T^ (? t E + VxVxE = S, (4.11) 

Here, S is a known term, representing the induced electromagnetic field of the excitation. - 
By vector analysis identities, (4.12) is equivalent with 

efiodttE + anodtV + V( V • E) - AE = S, (4.12) 

where A is the Laplace operator applied component- wise. 

Imaging: 

The initial pressure distribution can be obtained by acoustic reconstruction formulas. To 
calculate the conductivity a, we have now the quantitative aspect of thermo-acoustic to- 
mography: 

Problem 10 (qTAT). Let 

e^odttE + anadt'E + V(V • E) - AE = S 

with known S and constant e, (1q. 

Determine the conductivity a from knowledge of S and <r|E| 2 . 
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In [15], unique recovery of a from H = cr(x)|E(x)| 2 is shown for a class of bounded 
conductivity distributions (assuming that e = const., such that V-E = by the second 
eqn. in (2.1)). Also, the following scalar model is investigated 

EHodttu + a/iodtu — Au = 5, 

which is converted by the Fourier transform into 

Au + k 2 u 4- ikqu = S. 

Here, k = J— is the wave number, and a — —a\P^-. An exact reconstruction formula 
for q in this equation is developed in [8] (also assuming e is constant). 

4.2.3 Magneto-acoustic tomography with magnetic induction (MAT-MI): 

Magneto-acoustic tomography exploits the Lorentz force effect to create a pressure wave. 

Experiment: 

The measurement equipment consists of an excitation coil, a static magnet and ultrasound 
transducers. While the tissue is kept in the static magnetic field, electrical pulses in the 
excitation coil induce a eddy currents in the tissue. These, in turn, are affected by the 
Lorentz force due to the static field. The force causes the tissue to displace locally, and a 
pressure wave results. This wave reaches the boundary, where ultrasound transducers record 
the signal, either at the surface of the body or in a contact medium, such as water. 

Experimental setups with different arrangements of induction coils are documented in 
[44,109]. 

Mathematical modeling: 

We model the ultrasound wave motion with the wave equation with outer force f . Applying 
(2.2) with vanishing charge density p — to the static field B = Bo, we obtain the expression 
for the Lorentz force density f = J x Bo (see also [88]). Following the discussion in section 
2.4, the quantitative description of the acoustic wave is governed by 

^a„p(x, t) - Ap(x, t) = V • (J x B ) + ^d t H, 

where H quantifies the absorbed energy due to eddy currents. Notwithstanding, the energy 
absorption effect can be ignored in the experiments of [109]. There, magnetic induction is 
used for excitation of electrical currents, and this entails smaller field strengths. Considering 
the usual scale of conductivity in biological tissue, the pressure due to the thermo-acoustic 
effect can be neglected compared with pressure due to Lorentz force 3 (see also [109, p. 5179]). 
Therefore the model is 

\d tt p(x, t) - Ap(x, t) = V • (J x B ), (4.13) 

c z 

together with the initial and boundary conditions 

p(0) = 0, d tP (0) = 0. (4.14) 

We assume that the time- and space variables can be separated. The inhomogeneity f 
is therefore V • f(x)gr(t), where f is the force density —J x Bo = — cr(x)E(x) x Bo(x). Due 
to the different frequency ranges of the electrical and the acoustical pulse, we assume again 
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that g is a delta distribution. According to Duhamel's principle (2.36), the inhomogeneous 
problem (4.13) and (4.14) can be converted to a Cauchy problem for the homogeneous wave 
equation [7,9]: 

\dup-Ap = (pt,t) e fi x (0,oo) 

d tP = -V- (ctE x B ) t = 
p = Q t = 0. 

Imaging: 

The first step of the reconstruction of the source term V • (crE x Bo). In contrast to the 
techniques of impedance-acoustic tomography and thermo-acoustic tomography, we now 
look for the initial velocity of the pressure, which can be recovered as well. 

As quantitative problem, it remains to recover the distribution of the conductivity a over 
fi. For this purpose, one has to describe the electrical field E. For this, we use the eddy 
current model (2.27) discussed in section 2.3, adapted here for the pulsed excitation current 
(see also [70]). Therefore, the hybrid imaging problem is 

Problem 11 (Inverse problem MAT-MI). Let 



V • (crVu) = -<9 t A p • Vct in fi 
d n u = — <9t A p • n on 9fi 



u dS = 0. 

on 



The vector potential A p |n on M 3 is determined from the excitation current current Ji\p by 

Mo f Ji(y,t) , 3 



47r7 R 2 |x-y| 

Given are Bo and V • (a(dtA p + Vtt) x Bo), from which to determine a. 
4.2.4 Magneto-acoustic-electrical tomography (MAET): 

In MAET, pressure waves excite the tissue, and a voltage is produced due to the Lorentz 
force effect. 



Experiment: 

The measurement equipment in MAET consists of a permanent magnet, a waveform gener- 
ator, and of electrodes to measure the voltage at the boundary. The sample is put into the 
static magnetic field, and expopsed to different ultrasound waves. Due to frictional forces, 
the ions in the tissue move. As there is a magnetic field present, the electrical charges sepa- 
rate by the Lorentz force mechanism. A voltage results, which is detected on the surface of 
the probe. - The experiment is repeated with different arrangements of the static field. 

The physical effect underlying this technique has been studied in [76,104], further imaging 
experiments have been reported in [39]. 

Modeling: 

To model MAET, we follow the approach in [61]. We assume that the static magnetic field B 
in the experiments is homogeneous, that is, B(x) = B. The Lorentz force exerted upon the 
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moving particles gives rise to an impressed current density Jl which is approximately [61,76] 

Jz,(x,t) = cr(x)B x v(x,t), 

where v is the velocity field of the fluid. 

For each time t the electrical field inside the domain is described by an electrical 
potential, E(x, t) = — Vtt(x, t). In absence of any other free charges, the potentials can be 
described by the conditions V • J = and J • n = for the total current J = -aVii + Jl, 
that is 

V • (aVu) = - V ■ (o\B x v), x in fi 

, x (4-16) 
d n u = — (B x v) ■ n, x on dfl 

The voltage measured at the boundary gives knowledge of u(x, t)|go. One can also 
assume that only a weighted average of u is known, that is 

M(t)= I I{x)u{x,t)dS{x). (4.17) 
Jan 

Mathematical analysis: 

Using identities from vector analysis, it can be shown that the measurements in (4.17) are 

M(t) = B / J,xvd 3 x, (4.18) 
Jn 

with the current density J/(x) = — <tVwj(x) satisfying the electrostatic equation 

V • (erVu/) = 0, x in 

<JO n ui = i(xj, x on aS2. 

The inverse problem is 

Problem 12 (MAET). Lei satisfy (4.19). Determine a from knowledge of M(t) in 
(4.18), by possibly varying B, /, and v. 

Reconstructions for the case of general acoustic waves have been given in [61]. The 
procedure leads to knowledge of several vector valued current densities J/, which in the last 
step are converted to a. This is a problem of current density impedance imaging (CDII), 
discussed in section 4.1.1. Reconstructions for the case of focused ultrasound excitation are 
given in [7] . 

4.2.5 Acousto-electrical tomography (AET): 



An EIT experiment is performed and, in addition, the conductivity values are perturbed via 
ultrasound excitation in parallel. 

Experiment: 

For the measurements, one needs the typical EIT electrode setup and a device to produce 
ultrasound waves to excite the probe. The standard EIT measurement protocol is performed. 
For the perturbation of the conductivity, a set of various ultrasound waves is applied. The 
boundary voltages due to the perturbation are then compared with the standard voltages. 
The method was proposed in [110]. The underlying physical effect is treated in [66]. 
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Modeling: 

The standard electrostatic equation from (2.20) is used to model the electrical potential u, 

V ■ (txVu) = in D, 

a d n u = g on dQ. 

We now assume that the probe is excited by focused ultrasound beams. Therefore, the 
conductivity (which depends on the pressure), is perturbed in a domain D = z + SB, where 
B is the unit disc and z is the center of focus. The perturbed conductivity is denoted erf; 
the new potential u s is therefore described by 

V • (crfVuf) = in 

trf d n u s z = g on dtt. 

As measurement data, we have the difference of the measured voltages at the boundary, i.e. 
u d - u\ d Q. 

Mathematical analysis: 

It is shown in [5], that for S — >• 0, one has the asymptotic formula 

/ (u s -u)gdS=\Vu(z)\ 2 f cr(y) (??(y) ~ 1)2 rf 3 y + o(vol(£>)). 

Jen Jd v(y) + 1 

Here, rj is a known function in fi. 

From this formulation, one notices that the measurement data give knowledge of the 
power density at z, cr(z)| Vw(z)| 2 , up to a small error. 

Imaging: 

The inverse problem is to recover the conductivity er from the interior information er|Vu|. 
This problem is identical with problem 9 of impedance- acoustic tomography. 

Problem 13 (AET, conductivity determination). Let 

V ■ (erVit) = in 9. 

u\an = f ondfl 
a d n u = g on d£l. 

Given are f, g, as well as interior data H = cr|Vu| 2 . Determine a. 
Mathematical analysis: 

In the modeling of AET, we have assumed that cr z —a are approximately delta-distributions. 
But this restriction can be overcome and other more practical forms of ultrasound pertur- 
bation (such as plane waves, e.g.) can be used [57], [58, Section 2]. 
For this we consider 

M/)(y)= / fc(x,y)/(x)d 3 x 
Jn 

where L a maps an increment / = da of conductivity to the corresponding perturbation 
L a f in the boundary voltage. As perfect focusing of ultrasound waves correspond to 5- 
distributions /, and this would correspond to simple evaluation of the kernel k(x,y). It 
is more realistic that ultrasound waves of other, unfocused type, trigger the corresponding 
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Table 3: Electromagnetic variants of the Calderon problem in different frequency ranges. 



Technique 


equation in il 


boundary data 


remarks 


EIT 
eqn. (1.1) 


V -aVu = 


u 

ad n u 




IC-EIT 
section 4.1.3 


V • (crVu) = -iojA ■ V(J 


u 

a d n u = — iuiA. ■ n 


excitation A known 


MIT 
section 4.1.4 


V- (aVu) = -iuA-Va 


Jc n - B 

a d n u = — iuiA. ■ n 


excitation A known. 
B = B(ct,u,A), 
C is a coil surface 
outside f2. 


MD-EIT 
section 4.1.2 


r> r JxR 


B 


recover J, not a 



conductivity changes, proportional to the pressure. So, the intermediate step (synthetic 
focusing) corresponds to finding the values of the kernel k(x,y) by measurements of 

fc(x,y)w Q (x)d 3 x 

for a particular class of waves w a . Synthetic focusing for spherical waves, plane waves has 
been discussed in [57, Section 2] ; numerical calculations are performed in [58] . 

5 Summary 

The hybrid imaging techniques and mathematical models discussed above are summarized 
in the following tables 4,5. Those, consist of four columns, summarizing the imaging data, 
the overspecified data of the quantitative inverse problem, and the equations of the forward 
problems. In mathematical terms, in general, these are inverse problems with interior data. 
The table, however, does not summarize unique idcntifability and unique reconstructabil- 
ity of the quantitative imaging data. Hybrid models coupling magnetic resonance effects 
with current densities are listed in table 4. Some of these, like MREIT/CDII, rely on the 
low-frequency description V • (crVtt) = 0, others use radiofrequency currents, where the 
electrostatic approximation is not available and quasistatic approximations are used. Cou- 
pling hybrid techniques with ultrasound modalities are summarized in table 5. They also 
range from low-frequency or quasistatic to high-frequency ranges. 4 

In addition to coupled Physics (hybrid) techniques with interior information, we also 
summarize the standard EIT problem (the Calderon problem) together with techniques 
using electromagnetic coupling in several frequency ranges (see tables 3). The corresponding 
reconstruction cannot be separated into two imaging problems with interior data. 



Shortly, we also mention some hybrid imaging techniques using other physical parameters, which nev- 
ertheless are formally related to the models using the electrostatic regime, V • (ctVu) = 0: 

tn quantitative photoacoustic tomography (qPAT) in the diffusion regime, the diffusivity k is recon- 
structed, while the absorption /i in the elliptic equation V ■ (kVm) + fiu = is known. The 
imaging data are H = /il [17,24], 

In microwave tomography, the generalization of AET, one recovers the parameters a and q in 
V • (aVu) + k 2 q u = 0, with a\ Vu| 2 and q\u\ given [6, 9]. 

In pptoacoustic tomography one also recovers a in V ■ (crVu) = with <r| Vu| 2 given, where a is the 
optical absorption [16]. 
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Table 4: Coupled Physics models with interioir data, based on MRI. 



Technique 


interior data 


equation in 17 


boundary data 


remarks 


CDII 
section 4.1.1 


a Vu or 
a\Vu\ 


V- (dVu) = 


u 




MREIT 
section 4.1.1 


B z 


V- (aVu) = 


a d n u 


B z = F{aVu,a) 


CDII-MI 

section 4.1.5 


crE 


V • (aVu) = -iwA • Ver 


a d n u — iuA. ■ n 


excitation A 
known 

E = — iujA — Vu 


ICMREIT 
section 4.1.5 


B z 


V • (itVu) = - iujA ■ Va 


a d n u = iuA. ■ n 


excitation A 
known 

B z = F(aVu,a) 



Table 5: Coupled Physics models with interior data, based on Ultrasound. 



Technique 


interior data 


equation in 17 


boundary data 


remarks 


IAT 
section 4.2.1 


a\Vu\ 2 


V • (ctVu) = 


u 

crd n u 




qTAT 
section 4.2.2 


cr|Vu| 2 


£/io<9 tt E + anod t E+ 
V x V x E = S 


S\ F 


Fnn = 


MAT-MI 
section 4.2.3 


V • (o-E x B ) 


V • (o-Vu) = —d t A ■ Va 


ad n u = —d t A ■ n 


excitation A 
known, 

E = -d t A - Vu 


MAET 
section 4.2.4 


B • j crVu/ X v 


V- (ffVuj) = 


Ul 

ad n uj = I 


B • J Q aVu, x 

v = Ln IudS 


AET 
section 4.2.5 


a\Vu\ 2 


V • (ctVu) = 


u 

ad n u 
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